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Abstract 



The STAR data on the multiplicity correlations between narrow pseudorapidity 
bins in the pp and AuAu collisions are discussed. The PYTHIA 8.145 generator is 
used for the pp data, and a naive superposition model is presented for the AuAu data. 
It is shown that the PYTHIA generator with default parameter values describes the 
pp data reasonably well, whereas the superposition model fails to reproduce the 
centrality dependence seen in the data. Some possible reasons for this failure and a 
comparison with other models are presented. 
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1 Introduction 

The STAR experiment at the RHIC accelerator has measured the dependence of the mul- 
tiplicity correlations between the symmetric narrow pseudorapidity bins as a function of 
pseudorapidity distance for pp and AuAu data at 200 GeV in the nucleon-nucleon center- 
of-mass (CM) system [TJ. The data were taken separately for different centrality classes in 
AuAu collisions. The authors noted that the HIJING model [2J describes reasonably well 
the pp data, but fails completely for the central AuAu data, whereas the predictions of the 
Parton String Model [3] lie closer to the AuAu data, but fail to describe the dependence 
on pseudorapidity distance in the pp data. 

The AuAu data were discussed in the subsequent paper by Y.-L. Yan et al. [4j, who 
compared them with the PACIAE model [5 J based on the PYTHIA generator [Bj. The 
model overestimates systematically the correlation strength, and the difference increases 
strongly with decreasing centrality. 

The authors show also the predictions of three other models, which differ significantly 
from those of PACIAE model for most central events, but overestimate the correlations 
strength for semiperipheral events as strongly as PACIAE. 

In this paper we compare the STAR data and the models mentioned above with the 
new version of PYTHIA generator [7J for the pp data and the naive superposition model 
for the AuAu data. We have found that already with the default values of parameters 
the PYTHIA 8.145 generator describes reasonably well the pp data. On the other hand, 
the superposition model fails completely to describe the AuAu data. In particular, the 
centrality dependence of the correlation strength, which was weaker than in the data in 
the PACIAE model, is completely absent in our model. 

In the next section we describe the procedures used in the STAR paper for pp and 
AuAu data, which subsequently are used by us to analyze the events generated by 
PYTHIA and the superposition model. In Section 3. we show the results, and we discuss 
them in Section 4. Some conclusions and perspectives are presented in the last section. 

2 Definitions and procedures 

A standard quantity describing the strength of correlations between the phase space bins 
1 and 2 is the correlation coefficient b. 
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where n\ and 7i2 denote the multiplicities in bin 1 and 2, respectively. In the STAR data 
analysis the pseudorapidity bins of the width 0.2 (placed symmetrically in the CM frame) 
were used. The distance between the bin centers ranged from 0.2 to 1.8. 

It is well known that the correlation coefficients in hadron-hadron scattering increase 
with energy reflecting mainly the fact that the multiplicity distributions are much wider 
than Poissonian and their widths grow with energy. The STAR authors were mainly 
interested in the comparison of the correlation strength for different centrality classes in 
heavy ion collisions. These classes were defined by the range of multiplicity N of charged 
particles with p? > 0.15 GeV/c in a control bin of unit width, disjunctive with the bins 
for which b was measured (to be defined in Section 3). 

The b values calculated directly from the formula above reflect mainly the spread of 
the average multiplicities within each centrality class, which is quite large. The ranges of 



TV producing equal number of events in each class are approximately > 430, 320 — 430, 
230-320, 155-230 and 90-155 for the most central 0-10%, 10-20%, 20-30%, 30-40% 
and 40 — 50% events, respectively. The relation between n~l and A is approximately given 
by the ratio of the bin widths: nl = 0.2A, since dN/dn is nearly flat. 

Let us consider a distribution of a variable n depending on the external parameter A 
in such a way, that n = aN. It is well known that in such a case the dispersion may be 
separated into two parts: the first one averaging over this parameter, and the second one 
reflecting the spread of this parameter. The corresponding formula reads as follows: 

D 2 =<^2>-<n> 2 = < n^-n 2 > + <n 2 > - <n> 2 = < D 2 (N) > +a 2 D 2 N . (2) 

Here the bar denotes averaging over the distribution of n for given A, < ... > denotes 
averaging over N, and in the last term D 2 N denotes the dispersion of A. 

The authors of STAR performed a more involved procedure to separate the "dynam- 
ical" part of the correlation coefficient from the dominant "combinatorial" part. In each 
centrality class and for each distance between bins they calculated the averages nl, n 2 and 
n\ ■ n 2 for each value of A and fitted the dependence on iV within each class by the linear 
or quadratic formula: 

nl = di + bi-N, rf = c t + di-N + e 4 -A 2 , n x ■ n 2 = c 12 + d 12 ■ A + e 12 ■ N 2 . (3) 

Using the fitted values of the parameters they calculated these averages for iV = A 
and inserted them into the formula ([T|). A value of b was obtained for each centrality class 
and each value of pseudorapidity bin separation. 

One may note that the distributions of n^ for fixed A are narrow. We assume that 

D 2 (N) ~ 71(A). (4) 

It is easy to check that in such a case the procedure described above gives quite similar 
results as simple averaging of b over the distribution of A within each class. 

Indeed, let us assume that the formula (J4j) holds for both (identical) distributions of 
n\ and n 2 , and a similar formula holds for D 2 2 . Then the parameters of the fits ([3]) obey 
the following relations: 

a\ = «2, bi = b 2 , ci = c 2 , di = d 2 , e x = e 2 = e\ 2 = b\. (5) 

Averaging separately D 2 2 (N) and D 2 = D 2 = D\D 2 and inserting them into formula ([T]) 
we get in general 
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whereas from the STAR prescription we obtain 
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= _ D 2 12 (< A >) __ (ei 2 - b{) < A > 2 +(d 12 - 20^0 < A > +c 12 - a 2 

D 2 (< A >) (ei - b\) < A > 2 +(d t - 2oi6i) < A > + Cl - a\ U 

As we see, these two formulae are equivalent when the conditions fl5]) are fulfilled. 

For the pp data the STAR Collaboration measured the correlation coefficient b for each 
value of bin separation An and for each value of the multiplicity A in the pseudorapidity 
range of two units. Then the result was averaged over the distribution of A. This amounts 
approximately to using the formula (jH]). In fact, all three prescriptions are approximately 
equivalent. 



3 PYTHIA, superposition model and the data 

To understand the physical meaning of the data presented by the STAR Collaboration 
it is useful to perform the same procedures on the events generated by a Monte Carlo 
generator. For pp collisions there are many models describing high energy data. We use 
here the PYTHIA 8.145 with the default value of parameters. We analyze the generated 
events in the same way as described above for the STAR data. We checked that the 
PYTHIA results are practically the same when we use the non-diffractive (ND) and the 
non-single-diffractive (NSD) samples of events. 

In Fig.l the STAR data at Ecu = 200GeV and the model predictions for pp collisions 
are shown. The HIJING [2] and Parton String Model (PSM) [3] predictions are copied 
from the STAR paper pQ. The values of b depend rather weakly on the bin separation At] 
and the HIJING and PYTHIA models describe them reasonably well, although the A77 
dependence seems to be too strong in HIJING and too weak in PYTHIA. The PSM model 
predicts practically no Arj dependence, and seems to be incompatible with the data. 
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Figure 1: The correlation coefficients for the symmetric pseudorapidity bins of width 0.2 for the distance 
between the bin centers ranging from 0.2 to 1.8. Full circles in open squares, squares, triangles and full 
circles denote data and the predictions of PSM, HIJING and PYTHIA models, respectively. Model 
predictions are connected by lines to guide the eye. 



For AuAu collisions the situation is more involved. The centrality classes were defined 
by the multiplicity in a control pseudorapidity bin of the unit length, disjoint with the bins 
between which the correlation is measured. For the correlated bins (—1, —0.8), (0.8, 1) 
and (-0.8,-0.6), (0.6,0.8) the control bin is (-0.5,0.5); for the bins (-0.2,0), (0,0.2) 
and (—0.4,-0.2), (0.2,0.4) the control bin is made up of two parts: (—1,-0.5) and 
(0.5,1). For the bins (—0.6,-0.4), (0.4,0.6) the control bin consists of three parts: 
(—1,-0.75), (—0.25,0.25) and (0.75,1). All the generated samples of events were di- 
vided into the equally populated centrality classes, as described in the previous section. 

To generate a sample of AuAu events comparable with the data we used a naive super- 
position model. In this model a heavy ion collision with N part nucleon participants from 
each nucleus is presented as a straight superposition of N part pp collisions, each generated 
by the PYTHIA. To reproduce a sample of events belonging to the fixed centrality class 
we use the known experimental multiplicity distribution in the control bin P(N). The 
shape of the distribution of N part is assumed to be the same as that of P(N), and the 
average value < N part > is chosen to reproduce the value of < iV > in this class. 



In fact, the shape of P{N) is nearly linear in N for the large part of the distribution. 
Only for the most central part of the distribution (0 — 10%) it falls down much faster (ap- 
proximately Gaussian-like), and for the most peripheral events the maximum is sharper. 
This makes the generation of the P(N part ) quite easy. In Fig. 2 the predictions from such 
a model for five classes of the most central events are presented on the right hand side 
plot. On the left hand side the STAR data for the same centrality classes are presented 
and compared with the results of the PACIAE model [5j . 
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Figure 2: The correlation coefficients 6 as a function of pseudorapidity distance A77. On the left: the 
STAR data (full symbols) and PACIAE model (open symbols). Squares, circles, straight and inverted 
triangles and diamonds are used for the centrality classes 0-10%, 10-20%, 20-30%, 30-40% and 40-50%, 
respectively. On the right: our superposition model. The results for different centrality classes (marked 
as for the PACIAE model) are indistinguishable. 

We see that none of the models describes the data properly. For the most central 
class of events the disagreement of the PACIAE model with data is not very bad, but 
it grows quickly for the more peripheral classes. The Arj dependence of the correlation 
coefficient is rather weak and similar in the model and data. However, the dependence 
on the centrality is quite strong in the data and much weaker in the model. Our naive 
model shows practically no dependence on centrality. Thus the disagreement with data is 
significantly worse. The Arj dependence is similarly weak in the data and in the PACIAE 
model. 

One may add that the authors of the PACIAE model compared the data also with 
other models [5 J . The spread of their predictions is quite large and for most central events 
they bracket the data. However, the results for most peripheral data considered (the 
40% — 50% class) for all the models are much higher than the data, which exhibit in fact 
the negative values of the correlation coefficient b for the most distant bins. 

4 Discussion 



We have seen that the pp data analyzed separately for the given multiplicity N in the 
available pseudorapidity range and then averaged over the distribution of N are reasonably 



well described by different models. Both the value of the correlation coefficient 6 Ri 0.1 
and the weak dependence on the psudorapidity separation Ar] appear without any tuning. 

We should remember that the procedure used by STAR was designed to remove the 
influence of the global multiplicity fluctuations from the correlation effects. In fact, for 
the standard inclusive definition of the correlation coefficient b the dispersions in the 
numerator and denominator of formula ([1]) are dominated by these fluctuations. The 
value of b is then twice bigger than 0.1 and increases with the bin width. Thus a valid 
question is: what are the correlations measured by the STAR procedure? 

To clarify this point, let us consider an oversimplified model where there are no corre- 
lations for a fixed value of N. Then the distribution of the multiplicity n inside a narrow 
pseudorapidity bin of the width y is given by the binomial distribution 

P(n) = ( N \ n {\ - p) N ~ n ; p= <n> /N. (8) 

The generating function of this distribution is 

g(z)=^2^P(n) = (pz + l-p) N . (9) 

n 
Remember that the factorial moments of the distribution (jSJ) are given by 

n\ dg(z 



<n q > 
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It is easy to generalize these formulae for the case of two bins. The corresponding gener- 
ating function is then 

g{z llZ2 )= J2 z i 1 4 2 P(ni,n 2 ) = [p(z 1 + z 2 ) + l-2p} N . (11) 

The factorial moments for both bins can be calculated as in equation (TTUj) and the average 
product of multiplicities as 



^ _ d 2 g(z 1 ,z 2 ) 
< n x n 2 > = 



= p 2 N(N-l). (12) 

21 = 1, Z 2 = l 



dz\dz 2 
This results in the simple formulae for the dispersions 

D 2 l = Dl=p{l-p)N 1 D\ 2 = -p 2 N (13) 

and for the correlation coefficient 

Since this formula does not depend on N, the averaging of b over the distribution of N will 
give the same result. This will not change if we average the numerator and denominator 
in (|T%|) separately. 

We see that this simplified model disagrees strongly with data: the values of b are 
always negative, whereas in the data they are positive. Possible effects which introduce 
the positive correlations are: 



• the existence of resonances, probably responsible for the dependence on At] 

• the existence of jets and minijets, possible dijet correlations 

• the contribution of diffractive events (important for small N, decreasing for large 

N) 

• Bose-Einstein interference, which may be important for adjacent bins. 

Therefore we may conclude that the STAR procedure allows to measure the influence 
of these effects. They are significant, as they reverse the sign of the correlation coefficient. 
As indicated above, the realistic MC models describe them quite well. 

For the AuAu collisions none of the considered models describes the data properly. 
The centrality dependence of the correlation coefficient b seen in the data is significantly 
stronger than in the models. Moreover, our naive model predicts practically no centrality 
dependence at all and thus the disagreement with the data is even worse than for other 
models. 

Therefore one should ask what is the difference between our model and the other 
models. In our model it is assumed explicitly that the number of participating nucleons is 
the same in both colliding nuclei. It seems that this assumption suppresses the centrality 
dependence of the correlation coefficient. 

This observation suggests that the data of STAR Collaboration, when interpreted 
within the superposition picture, measure mainly the fluctuations in the number of par- 
ticipating nucleons N part in both nuclei. For the most central class of events these fluctua- 
tions seem to be strongly correlated. With decreasing centrality this correlation becomes 
weaker both in the data and in the models. 

One can understand this effect qualitatively in the framework of a "hidden asymmetry" 
picture |8J, in which the dominant configuration is strongly asymmetric. There are few 
events in which the number of participants in both colliding nuclei is the same or similar. 
The correlation coefficient b is positive and large for most central events, where the number 
of participants in both nuclei is nearly maximal (and thus approximately the same). The 
centrality is defined by the multiplicity in a symmetric control bin close to r\ — 0. When 
this multiplicity decreases, the asymmetric configurations start to dominate, which results 
in decreasing values of b. This is illustrated in Fig.3, where our toy model predictions 
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Figure 3: The correlation coefficients b as a function of pseudorapidity distance A 77 for the 40 — 50% 
class. Asterisks, triangles and circles are for the original model, triangular and flat distributions of Np, 
respectively. 



for the 40 — 50% class are compared with its two modified versions. In these versions 
the particles from each p — p event generated by PYTHIA are divided randomly into the 
fragmentation products of two "wounded" nucleons. One of such "half-events" populates 
mainly the forward CM hemisphere and the other one backward one. The AuAu event 
is constructed from the (generally unequal) numbers of such half-events Np and Nb- For 
a given sum of these numbers we tried the distribution of Np to be flat or to have a 
triangular shape. 

We see that it is quite easy to explain the decrease of b with the pseudorapidity 
distance, even to the values much more negative than in the data. However, our model 
fails to explain the strong increase of b for most central events and the adjacent bins. 

This dependence may be reproduced qualitatively in PACIAE and other similar mod- 
els, because they are not "pure" superposition models. The partons from initial nucleons 
interact forming the parton cascades which then hadronize. With increasing centrality 
the number of parton cascades increases. Apparently, the corresponding increase of the 
numerator in formula (pQ) is faster than the corresponding increase of the denominator. 

However, the decrease for non-central events is underestimated in all the models. The 
negative values of b in the data for the 40 — 50% class, which correspond to the fact 
that for this class of events the positive fluctuation of N part in one nucleus is correlated 
with a negative fluctuation in the second nucleus, are not reproduced. This conforms 
with the idea of "hidden asymmetry". If asymmetric configurations are prevalent for 
semiperipheral events, the anticorrelation is natural. 

5 Conclusions and outlook 

The failure of superposition models is usually explained by the collective effects. An 
example of such an effect is the quark-gluon plasma formation, which should be most 
visible in the central events. However, the data for the most central 10% of events are 
quite well described by the models where no plasma is formed |4]. As noted above, the 
biggest discrepancy is observed for semi-peripheral events. Other collective effects, as the 
hydrodynamical flow in the hadronic phase, are strongest for this class of events. On the 
other hand, it is not clear how such effects may be responsible for the discrepancy between 
the models and data. It seems more likely that the explanation is related to the "hidden 
asymmetry" of the interacting nucleons. 

The negative correlation for the semiperipheral collisions are not the only unexpected 
feature of the STAR data. It was noted [9J that the large values of b for the most 
central class of events suggest that the correlation between the distant bins "B" and "F" 
is stronger than the correlations between these bins and the central control bin. This 
contradicts the generally accepted feature of the uniform decrease of correlations with 
increasing distance between bins. Therefore we conclude that the centrality dependence 
of b reported by STAR seems to be surprisingly strong. It would be interesting to measure 
this dependence collecting further data. 
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